Triangle Anomalies, Thermodynamics, and Hydrodynamics 

Kristan Jensen 1 

1 Department of Physics, University of Victoria, Victoria, BC V8W 3P6, Canada 

(Dated: February 27, 2013) 

We consider 3+1-dimensional fluids with U(l) 3 anomalies. We use Ward identities to constrain 
low-momentum Euclidean correlation functions and obtain differential equations that relate two 
and three-point functions. The solution to those equations yields, among other things, the chiral 
magnetic conductivity. We then compute zero-frequency functions in hydrodynamics and show that 
the consistency of the hydrodynamic theory also fixes the anomaly-induced conductivities. 
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Introduction. — The chiral magnetic effect (CME) [l|-|3| 
and chiral vortical effect (CVE) [J, |5[ represent remark- 
able implications of anomalies in quantum field theory for 
macroscopic transport. In a fluid with a U(l) 3 anomaly, 
there will be currents directed along a magnetic field B^ 
or the vorticity w^ . The former is the CME and the latter 
the CVE. Both effects probe the violation of parity and 
charge conjugation and so may be measured by studying 
spatial and charge asymmetries [5| in off-axis heavy ion 
collisions at RHIC or the LHC. 

It has also been shown that the hydrodynamic descrip- 
tion of a relativistic fluid with anomalies must be mod- 
ified [4|] from its textbook treatment [6J. There are ad- 
ditional transport coefficients, describing the response of 
the currents to magnetic fields and vorticity. These co- 
efficients are fixed by demanding a local version of the 
second law of thermodynamics, namely the existence of 
an entropy current whose divergence is positive semi- 
definite |6|. Despite the absence of a clear understand- 
ing of this constraint in field theory, the result for the 
anomaly-induced transport matches results at weak [2j 
and strong coupling [3, [8| . This yields a hydrodynamic 
description of the CME and CVE. 

We endeavor to reproduce the constraints on anomaly- 
induced transport without recourse to the entropy cur- 
rent. We will do so by employing properties of equi- 
librium quantum field theory. In particular, we will 
study theories which have finite static correlation length, 
A. This has the practical implication that real-space, 
zero-frequency correlation functions fall off exponentially 
at long distance. The Fourier transformed functions 
are then analytic at zero frequency and small momenta 
k\ -C 1. Together with some Ward identites we relate 
the U(l) 3 anomaly coefficient to transport. 

The resulting calculation makes a few points clear. 
First, only the small B^ and small w^ parts of the CME 
and CVE are fixed by the anomaly. Second, anomaly- 
constrained transport follows from a covariant and gauge- 
invariant description of thermodynamics. Finally, we see 
that, among other things, the entropy current enforces 
properties of the equilibrium theory (see ([5]) and (|SJ|) that 
are not manifest in the hydrodynamic description. 

Note: We explore a number of questions related to 
equilibrium thermodynamics and the role of the entropy 



current in a companion paper [9|. While this work was 
in progress we were also made aware of 10] which has 
overlap with the content of this Letter. 

Correlation functions atT,fj, ^ 0. — In real-time finite 
temperature field theory there arc different definitions for 
correlation functions. We employ the closed-time-path 
(CTP) formalism, in which time is extended to a closed 
contour which extends from t\ : (— oo, oo), and then dou- 
bles back as t% : (oo,— oo). See [11| for a review. The 
only ingredient we need is that there is a CTP generating 
functional Wctp which depends on two sets of sources: 
</,. = (Ji + J<i)j1, J a = J\ — J2, where the Ji are inde- 
pendent sources introduced on each segment of the time 
contour. The fully retarded functions, which are those 
computed in hydrodynamics [12], are obtained by vary- 
ing Wctp with respect to the r sources and one a source. 
The r source couples to a-type operators, whereas the a 
source couples to r-type operators. Insertions of U{1) 
currents, labeled by the index A, and stress tensor densi- 
ties come from varying Wctp with respect to background 
fields as 
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The currents and stress tensor are related to the densities 
by J^ A = y/^gJ^' A ,T^ v = y/=gT» u . In this work we 
consider ra and raa functions in momentum space. We 
notate these as 



G I /(q) = (0 I r (q)0 J a (-q)), 
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( qi ,q 2 ) = (OUqJO'MO* {-qi - q 2 )), (3) 



where coordinate space functions are related to their mo- 
mentum space cousins by 

f(xi,..,x n ) = [d d qi ..d d qn e^^ + -i^f( qi ,..,q n ). 



_ (4) 
We further specialize to zero frequency ra..a functions, 
i.e. we take qi = (0,k;). These (and all other CTP func- 
tions) are proportional to the corresponding Euclidean 
time-ordered functions obtained by variation of the Eu- 
clidean generating functional We UM- Because of this 
we henceforth neglect the r and a indices. Additionally, 



since bosonic derivatives of We commute, the CTP func- 
tions of bosonic fields satisfy 
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The diffeomorphism invariance and anomalous varia- 
tion of Wctp leads to the Ward identities 14[ , 
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where C ABC is the symmetric anomaly coefficient and 
£ oi23 = i , (jgj an d (|7J) encode local gauge and diffeomor- 
phism covariance. There are additional Ward identities 
for the Euclidean theory at nonzero temperature due to 
the global topology of R 3 x S 1 15]. Constant A Q and 
hoo may be gauged away at the price of redefining the 
temperature and chemical potentials. These are shifted 
as 
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variation of ([6]). Setting the two equal gives 
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By applying the same method to any Gr with a J^ 
insertion we thereby find 
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Imposing (J7J fixes the remaining S's. For example, 
^^'^(ki.k,) = -e tt, (* 1 ) fc (fc)i£!l£ fl (13) 

J 2 (^r^°(k 2 )) 
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= i(ki+k 2 ) k G R ' % (-ki), 



The /i" 4 depend on goo because they are defined through 
the Wilson line of A A around the time circle. 

At small momenta |k,|A <C 1, the correlation functions 
of the current and stress tensor are heavily constrained. 
The two-point functions with parity-violating terms to 
0{k) may be parametrized as 
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for some functions a m . By ^ a AB — a BA . The second 
comes from the Bardeen-Zumino poly- 



term in G 
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nomial, which encodes the anomalous dependence of the 
jii>A on g au ge fields 16j. The three-point functions with 
parity-violating terms to 0(k) are similarly constrained. 
Demanding that the three-point functions are consistent 
with ([3]) we find that they take the form 
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Imposing (J6|) fixes half of the E's. For example, we 
may compute i(k\)iG R ' 3 ' directly from ((TU)) or by 



which by © and ([5]) gives 
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Applying this method to the other three-point functions 
in (fTU|) yields the remaining S's 
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By the discussion around © we may evaluate the two- 
point functions © in a background with constant A A and 
goo- To 0(k) we find 
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where the prime indicates that a quantity is evaluated 
at temperature T" and chemical potentials /i' A (|HJ|, and 
the subscript S* that the correlator is evaluated in the 
presence of background fields. 

Differentiating (fT6|) with respect to A A and 500 leads 
to three-point functions with zero momentum insertions 
of J°' A and T 00 . Comparing these functions with the 
three-point functions (ITU1) and using (fT2")l , (|T4")) , and (|15l) , 
the two agree only if the six equations 
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are satisfied. These equations uniquely fix the a m up to 
integration constants. We have 
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We may also consider the behavior of the integration 
constants f m under CPT. By a hydrodynamic argument 
employed in [l7| , we find that the f AB and / 3 are CPT- 
violating, while the f A are CPT-preserving. We can also 
directly establish this result by studying the transforma- 
tion of the two-point functions (j9]) under CPT. 

Hydrodynamics with sources. — It is instructive to re- 
produce (|19|) from hydrodynamics. In this section we 
begin with equilibria at constant T, p , and vanish- 
ing sources. We take the fluid rest frame to be u^ — 
v^/y— v 2 , with v p a constant timelike vector. In these 
states the stress tensor and current are 



(T"") = eu*V + PA"", (J»- A ) = p 
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with P and e the pressure and energy density, obeying 

dP = sdT + p A dp A , e + P = sT + p, A p A , (21) 

and A pu = g pu + u^u v is a projector satisfying A p "u v = 
0, A 2 = A. In hydrodynamics we study long- wavelength 
fluctuations around equilibrium states. Those fluctua- 
tions may be described by promoting T, p , and the u p 
to spacetime fields (the hydrodynamic variables) and ex- 
panding the stress tensor and current in gradients thereof 
- this is the derivative expansion! [18| . We also turn on 
slowly varying background gauge fields A A and metric 
perturbations g — rj + h. We take the sources to be O(l), 
so that the field strengths F and connection coefficients T 
are 0(d) in the derivative expansion. This is the scaling 
required to study the response of a fluid to sources. 
In more general states we have the decomposition 

(T"") = £u^u v + VA» V + q p u v + q v u^ + t" u , 

ftABC 

U^ A ) = N A u» + v^ A + —^e^ p,y A B FZ, (22) 
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(f2"2"j) has some redundancy. Fixing (T) and (J A ), we may 
redefine the hydrodynamic variables by 0(d) quantities 
to impose a choice of hydrodynamic frame. In the rest 
of this work, we perturb the equilibrium state (|20|) by 
sources which are static in the fluid rest frame, supposing 
that the perturbed fluid remains in a time-independent 
equilibrium. In such a state the one-point functions (|2"2")l 
and hydrodynamic variables will be local functions of 
sources - the non-locality is on order of the size of the 





1 


2 


vectors (vf ) 


Ei 


a M 


psuedovectors (vf) 


J7== Q ^ P °^F A 


1 ^ vpa UvdpUa 



TABLE I. The independent first-order tensors. 



static screening lengths, which vanish in the hydrody- 
namic approximation. We then perform a change of 
frame such that T, p , and v p take on their form for 
equilibria with constant background fields, 
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where T cq and v£ are the temperature and velocity field 
of the source-free equilibrium state. This is the covariant 
version of (JSJ). 

It remains to express £ , V, A/", q^, v^ , and t^ v in terms 
of the sources. These are the constitutive relations. To 
proceed we compute derivatives of T, p , and the fluid 
velocity. We find 



dfj,T = -Tcifj,, 



d fi p A = -p A a fi + E A , (24a) 
V M u„ = -u^a v + W/M/ , E A = F A „u u , (24b) 

are identically satisfied with 
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a" = u u V„u», u» v = (V p zv - V ff u p ). (25) 

The independent tensors with one derivative are listed in 
Table HI The pseudovectors are v p ' = B^' A , the mag- 
netic field, and v^ = w^ the vorticity of the fluid. To 
0(d) we then write 



£ = e, V = P, Af = p, r"" = 0, 



(26) 



where the 7's, 7's, <5's, and <5's are functions of T and p A . 
We continue by treating the Ward identities ([6]) and ([7]) 
as equations of motion. Ordinarily, we solve for the hy- 
drodynamic variables in the presence of external fields. 
In this instance, the conservation equations leads to dif- 
ferential equations for the coefficients in the constitutive 
relations. We note that the gradients of the first-order 
tensors satisfy some simple relations 
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where R^ v is the Ricci tensor. Applying (|24f and (j27|) 
to the conservation equations leads to a sum of coeffi- 
cients, each of which multiplies an independent second- 
order tensor. Each such coefficient must vanish, which 
leads to a number of equations, 
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When (HH1) holds, the Ward identities © and are ex- 
actly solved for the first-order constitutive relations (f2l)| . 
We compute the ra..a functions by varying one-point 
functions (which we view as expectation values of r oper- 
ators) with respect to external fields |l2j (which we view 
as a- type sources). Those variations are easy to perform 
in this frame. We obtain zero-frequency n-point fuctions 
by directly varying (f2"2"]l and (|26[) . For instance, we find 
the following two-point functions to 0(k), 



Zero-frequency correlation functions encode the ther- 
modynamic response of a fluid. In this instance, a fluid 
may be subjected to a magnetic field B^ or vorticity w^ 
and remain in equilibrium. When parity is not a sym- 
metry, there may be charge and momentum currents di- 
rected along B^ and u> M . This calculation is telling us 
that the O(B) and 0(w) response of those currents is 
fixed by a consistent description of thermodynamics in 
the presence of background fields. 

One significant question remains: what is the role of 
the CPT-preserving integration constant f% (fl9f ? At 
weak [2JJ and strong 22| coupling, it has been found to 
be proportional to the mixed anomaly coefficient. Is that 
result fixed by Ward identities? The hydrodynamic cal- 
culation sheds some light on this question. By exactly 
solving the conservation equations to O(0 2 ) and impos- 
ing ([5]) at the outset, we have implicitly enforced the 
Ward identities for gauge/diffeomorphism invariance on 
the O(k) parts of all zero- frequency n-point functions of 
T and 3 '■ At the end of the day, f% is left unfixed by 
the zero- frequency conditions ([5|) and ((HJ. Three logical 
possibilities remain: (i.) f£ is fixed by another zero- 
frequency condition, (ii.) f% is fixed by a finite-frequency 
property like the Onsager relations, or (iii.) f% is unfixed. 
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precisely reproduce those (fTT)) 



so that the constraints ( 
and |TS]| for the cr's. 

Discussion. — In this letter we have used Ward iden- 
tities to constrain zero-frequency correlation functions. 
For a 3 + 1-dimensional theory with U(l) 3 anomalies, 
the O(k) parts of three-point functions of the stress ten- 
sor and currents are determined by conservation in terms 
of the anomaly coefficients and 0(k) parts of two-point 
functions. © leads to differential equations that relate 
two-point functions to three-point functions. As a re- 
sult the 0(k) terms in two and three-point functions of 
T and J are uniquely fixed up to integration constants. 
The result matches the hydrodynamic [J, [l_9j and holo- 
graphic [20( calculation up to additional CPT-violating 
integration constants f AB (fT9"l) . 



[1] D. Kharzeev and A. Zhitnitsky, |Nucl.Phys." A797, 67 

(2007) 
[2] K. Fukushima, D. E. Kharzeev, and H. J. Warringa, 

|Phys.Rev.| D78, 074033 (2008) 

[3] D. E. Kharzeev and H. J. Warringa, Phys.Rev. D80, 

034028 (2009), 10 pages, 4 figures 
[4] D. T. Son and P. Surowka, |Phys.Rev.LettJ, 103, 191601 

(2009) 

[5] D. E. Kharzeev and D. T. Son, |Phys.Rev.L"etE1 106, 

062301 (2011) 
[6] L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Perg- 

amon Press, 1987) 
[7] N. Banerjee et a UJHEPI Ol. 094 (2011) 
[8] J. Erdmenger, M. Haack, M. Kaminski, and A. Yarom, 

IJHEPI 01. 055 (2009) 
[9] K. Jensen, M. Kaminski, P. Kovtun, R. Meyer, A. Rtz, 
and A. Yarom(2012), larXiv:1 203.3556 
[10] N. Banerjee, J. Bhattacharya, S. Bhattacharyya, S. Jain, 

S. Minwalla, and T. Sharma(2012), arXiv: 1203.3544 
[11] E. Wang and U. W. Heinz, Phys. Rev. D66, 025008 

(2002) 

[12] G. D. Moore and K. A. Sohrabi, |Phys. Rev. Lett.| 106, 

122302 (2011) 

[13] T. Evans, |NucT.Phys.|B3 74, 340 (1992)B496, 486 (1997) 
[14] C. P. Herzog, |J.PJrJslA] A42, 343001 (2009) 



[15] K. Jensen, M. Kaminski, P. Kovtun, R. Meyer, A. Ritz, 

and A. Yarom(2011), arXiv:1112.4498 
[16] W. A. Bardeen and B. Zumino, |Nucl.Phys. B244, 421 

(1984) 
[17] J. Bhattacharya, S. Bhattacharyya, S. Minwalla, and 

A. Yaromf2011). larXiv:1105. 37331 
[18] S. Bhattacharyya, V. E. Hubeny, S. Minwalla, and 

M. Rangamani, JHEP 0802, 045 (2008) 



[19] Y. Neiman and Y. Oz, JHEP 1103, 023 (2011) 

[20] I. Amado, K. Landsteiner, and F. Pena-Benitez, JHEP 
1105, 081 (2011) 

[21] K. Landsteiner, E. Megias, and F. Pena-Benitez, 
|Phys.Rev.LettJ, 107, 021601 (2011) 

[22] K. Landsteiner, E. Megias, L. Melgar, and F. Pena- 
Benitez, JHEP 1109, 121 (2011) 



